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Abstract 

In this paper we prove a viability result for multidimensional, time dependent, 
stochastic differential equations driven by fractional Brownian motion with Hurst 
parameter^ < H < 1, using pathwise approach. The sufficient condition is also 
an alternative global existence result for the fractional differential equations with 
restrictions on the state. 
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1 Introduction 

Let B = {Bt,t > 0} be a fractional Brownian motion (fBm) of Hurst parameter 
H G (0, 1). That is, i? is a centered Gaussian process with the covariance function (see 

Rnis, t) = E (BA) = ^ (t^^ + - |t - . (1) 
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Notice that if if = i, the process i? is a standard Brownian motion, but if if 7^ |, it 
does not have independent increments. It clearly follows that K\Bt — Bgl"^ = |t — sp^. 
As a consequence, the process B has a— Holder continuous paths for all a G (0, H). 

The definition of stochastic integrals with respect to the fractional Brownian motion 
has been investigated by several authors. 

Essentially two different types of integrals can be defined: 

• The divergence integral (or Skorohod integral) with respect to fBm is defined as 
the adjoint of the derivative operator in the framework of the Malliavin calculus. 
This approach was introduced by Decreusefond and Ustiinel [7] and developed 
by Carmona and Coutin [6j, Duncan, Hu and Pasik-Duncan [8], Alos, Mazet and 
Nualart pLj, Hu and 0ksendal [10], among others. This stochastic integral can be 
expressed as the limit of Riemann sums defined using Wick products and satisfies 
the zero mean property 

• The pathwise Riemann-Stieltjes integral J^UsdB^ which exists if the stochastic 
process (wOtg[oT] continuous paths of order a > 1 — if, is a consequence 
of the result of Young [17]. Zahle has defined in [18] and [19] this integral for 
processes with paths in a fractional Sobolev type space. In this paper we will 
follow this last approach. 

The aim of this paper is to state necessary and sufficient conditions that guarantee 
that the solution of a given (forward) stochastic differential equation driven by the 
fractional Brownian motion B with Hurst parameter 1/2 < if < 1 (in short: f-SDE), 
P - a.s. u efl 

/sVt i>sVt 
b{r,X'^'^)dr + a{r,Xl'^)dBf, se[t,T], (2) 

{t,x) e [0,T] X evolves in a prescribed set K i.e., under which it holds that for all 
t e [0,T] and for x E K : 

Xl'^'eK a.s.ujen,^se[t,T]. 

where 

• B = [B^)i^^^ , B^, i = l,k, are independent fractional Brownian motions with 
Hurst parameter if, ^ < if < 1., and the integral with respect to i? is a pathwise 
Riemann-Stieltjes integral; 

• Xq = (Xq)^^-^ is a d - dimensional random variable defined in a complete proba- 
bility space (f2,JF, P); 

• 6 (cj, -, ■) : [0, T] X M^^M^ and a (cj, -, ■) : [0, T] x R'^'"' are continuous func- 
tions. 
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A general result on the existence and uniqueness of the solution for multidimen- 
sional, time dependent, stochastic differential equations driven by a fractional Brown- 
ian motion with Hurst parameter if > 1/2 has been given by Nualart, Ra§canu in [IQ\ 
using a techniques of the classical fractional calculus. 

The viability - property, has been extensively studied for deterministic differential 
equations and inclusions, starting with Nagumo's pioneer work in 1943 ( see [2] for 
references) . To our knowledge the first work that gives a characterization of the viability 
property in a stochastic framework was written by Aubin and Da Prato [3] in 1990 (see 
also [3], [9]). The key point of their work consists in defining a suitable Bouligand's 
stochastic tangent cone which generalizes the cone used in the study of the viability 
property for deterministic systems, (for other points of view see also [H], [13] and 

Another approach has been developed by Buckdahn, Quincampoix, Rainer, Ra§canu 
in [5j which present a unified approach for the study of the viability property of SDE, 
BSDE and PDE by relating the distance to the constraint to some suitable PDE. More 
precisely, in the case of an SDE (respectively BSDE) they show that the viability prop- 
erty is equivalent to the fact that the square of the distance function is a viscosity 
super solution (respectively subsolution) of the PDE. 

In this paper we will prove a Nagumo type Theorem on viability proprieties of close 
bounded subsets with respect to a stochastic differential equation driven by fractional 
Brownian motion, following an approach inspired by the work of Nualart and Ra§canu 

The organization of the paper is as follows. In Section 2 we recall some classical def- 
initions and consider the assumptions on the coefficients supposed to hold. In Section 
3 we state our main result. Section 4 contains the deterministic more general result 
and Section 5 contains the proof of the result in the stochastic case. 

2 Preliminaries 

2.1 Generalized Stieltjes integral 

Let d,k E N*. Given a matrix A = (a*'"')rfx/c ^^'^ ^ vector y = (i/*)(^xi denote 

1^1^ = Eij and = Y.i Iv'f ■ 

Let t e [0, T] be fixed. Denote by iy"'°°(t, T; M'^), < a< 1, the space of continuous 
functions / : [t, T] R"' such that 

ll/IU;[,T] := sup (\m\+ H-ll^l^I^dr) < oo. 
s€[t,T] V Jt [s-r) ) 

A equivalent norm can be defined by 

iii/iiL,M.,T, =' -.p (i/(»)i + r^/^'-ffi rfr) 
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for any A > 0. 

For any < < 1, denote by ([t, T] ; R*^) the space of /x- Holder continuous functions 
/ : [t, T] R'^, equipped with the norm 

, I/W-/MI 

t<s<T-<T 



f.;[t,T] = ll-^lloo;[t,r]+_SUp_^— -^<00, 



where ||/||oo-rtTi ^^P We have, for all < £ < o; 

se[t,T] 

C"+^([i, T] ; W^) C W'^^it, T; W^) C C°^-"([i, T] ; M'^). 

with continuous embeddings. 

Fix a parameter < a < |. Denote by W^~°''°°{t,T;W^) the space of continuous 
functions g : [t,T] such that 

I mi a- f \g{s)-g{r)\ , \g{y) - g(r)\ ^ \ ^ 

Clearly, 

C^-"+^([t, T] ; M'^) C T; R'^) C C^-°([t, T] ; R'^), > 0. 

Denoting 

Aa(^;[^,r])'=^— 1— sup |(i^lZ"y.-)(r)|, 

i (^i — t<r<s<T 

where 

r(Q;) = p°'''^e~Pdp is the Euler function and 

{D.^ (-■) = T(J^r lH^Tj^ ^ - / (y- rf-'V 

we have ^ 

Note that 

A„(^;[t,T])<A«(^;[0,T]) (=U^ig) 



Let T; R'') the space of measurable functions / on [t,T] such that 

a,l;[t,T] - I ~ ~ ' ~ r-T^^H ds < OO 



Clearly T; R<^) C W^^'^t, T; R<^) and \\fL,^.,^,,r] < (^^ + &) II/IL,oo;[,t] • 



Denote 
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Definition 2 Let < a < If f e VT^'^t, T; M'^^'^) and g G T; M'^), then 

defining 

f (r) dg (r) t/ (-1)" (r) (Dlz'^g,^) (r) dr. (4) 



t/ie integral fdg exists for all s E [t, T] and 

f f (r) dg (r) < sup | {Dlz'^g,,) (r) | f \ ( {s) \ ds 

Jt t<r<s<T Jt 



(5) 



< A.(^;[t,T]) 11/11 _ 



[t,T\ ■ 



We give two continuity properties of this integral used in the paper. Let t < s < 
T < T. We denote 



Gs,r {f) = J f (r) dg (r) = Gt,r (/) - Gt,s (/) 
Proposition 3 Let < a < ^. 

(a) Gt,T ■■ W'^it, T; R'^) ^^-"([t, T] ; R'^) is a Unear continuous map and 

\\Gt,- (/)lli-«;[,T] < <T A.(^7; [t,T]) ||/|L,^^[,^] (6) 

where A)^ j, is a positive constant depending only on a and T; A^J^^^ < 4 + 3T. 

(6) Gt,T '■ W°''°°(t,T;R'^) W°''°°(t,T;R'^) is a linear continuous map and for all 

A > 1 

where A^^j, is a positive constant depending only on a andT; A^^j, = - — — ^ — h . 
Proof, (a) Let t < t < s < T. From the definition of the integral we have 

and therefore 

(/)! < A^{g; [t, T]) (2 + T") (. - r)^^° H/L,,.,!,^] • (9) 

Hence the inequality ([6]) follows with A^^^^ = ^ + T + 2 + T" < 4 + 3T. 
(6) We have 
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\GtAf)-GtAf)[^^ 



f \ l+a 

[s — r 



t Jt 

S P0 



Since 



= (s - 0)-'° / (1 + u-'^du 
Jo 

(pi POO 



and fort<r<M<6'<s 



-2« 



/ (s — r)^"~ dr < — {s — m)^" < — (s — t) 
Ji a a 



-2a 

1 



then it follows 

\Gt,s{f)\+ / z rr+E, dr 

Jt [s - r) 

<K (^?;[t,T]) ( ^ 

\ (1 — a) a 



X ^' ((. - + (9 - r") (i/(»)i + /"^f=^* ) 

The inequality ([7]) clearly follows since 
J^'e-^(^-^) [(s - ^)-'" + - t)-°] rf^ 

1 4 



- Xl-2a I -2a 
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2.2 Assumptions and notations 

Consider the equation on R'' 



„5 m „g 

Xi = Xi,+ / b'ir, Xr)dr + V / o''^ (r, X,) dBl s E [0, T] 
Jo Jo 



(10) 



i = l,...,d, where the processes , j = l,...,k are independent fractional Brownian 
motions with Hurst parameter H defined in a complete probability space JF, P), 
Xo is a. d - dimensional random variable, and the coefficients a*'-', 6* : x [0, T] x M'^— >M 
are measurable functions. Setting a = (c"*'-')^^^ ' ^ ~ (^*)dxi ' -^s ~ i^s)kxi 
Xs = {XI) then we can write the equation ffTOl) in a simpler form 

PS ns 

Xs = Xo+ / b{r, Xr)dr + a {r, X,.) dBr, s G [0, T] . 
Jo io 

Let us consider the following assumptions on the coefficients, which are supposed 
to hold for P - almost all uj E Q. The constants Mq, M^j, Lq, may depend on u. 

(Hi) a{t, x) is differentiable in x, and there exist some constants (3,6, < (3,6 < 1, 
and for every R>0 there exists Mr > such that the following properties hold 
for all te[0,T],¥- a.s. cu e n : 



{Ha): 



II 



i) \a{t,x)~a{s,y)\<M^{\t~s\^+\x~y\), Vx, y G 

\V^a{t,y)-V^a{s,z) \ < MR{\t - s\f + \y - z\^) , ^\y\,\z\ < R, 



where Vxa{t,x) = (V^cr*'^ {'t,x))i=T4, j=Tji 

d d k 

|V.a(t,a;)|' = ^5^5^|a.y'^-(t,x) 



£=1 i=l j=l 

Remark that for all a; G M'^ 

\a{t,x)\ < |(T (0,0)1+ Mo (Itl'^ + lxl) <Mo,t(1 + |x|) 

where Mo,t = W (0, 0)| + Mq + MqT. 
Let 

ao = mm < -, p. 



2"^' 1 + 5 

With respect to the coefficient b we assume 

(H2) There exist yU G (1 — ao, 1] and for every R> there exists Lpi > such that the 
following properties hold for all t G [0, T] , P — a.s. u E fl : 



{Hb 



i) \b{r,x) -b{s,y)\ < LR{\r - sl^" + \x - y\) , V |a;| , |?/| < i?, 
ii) \b{t,x)\ < Lo{l + \x\), ^x eR'^. 
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From the work of D. Nualart and A. Ra§canu [16] we deduce that under the assump- 
tions (Hi) and {H2) with /3 > 1- H and 6 > - 1 and for every fixed (t, ^ G [0, T] x 
L° (n, J^, P ; R"^)) , the SDE has a unique solution G L° {n, P ; ^^"'""(t, T; R"^)) , 
for all a G (1 — -ff, ao) • Moreover, for P - almost all G 

X(.;,-)=(X^ gC^-" (0,T;R'^). 

In this paper we develop an other proof and we obtain a stronger existence result: if 
the starting point ^ belongs to a closed set K then there exists a unique solution 
X*'^ and the solution evolves in K. 



2.3 Basic estimates 



Fix a parameter < a < 2 ^ ^" ^^iven two functions / G T; R"') and g G 



Vr^-"'°°(t,r;M^) we denote 



Gg(/)= / a(r,/(r))ci^7(r) 



From Proposition [3] we infer 

Corollary 4 Lei the assumption (Hi) be satisfied and < a < | A /3. Let f G 

Vr"'°° (t, T; M^) an(i g G T; R^). Then G^^^ (/) G C^"" (t, T; R'^) an(i /or all 
A > 1 



l-a;[t,T] 



11] 



(/) 



a,A;[t,T] 



A 



l-2a 



a,A;[t,T] / ' 



where C^"^^^ and Cq*^^^ are constants which only depend on Mq^T; Mq, T, a, (5. 
Proof. From Proposition [3] we have 



Gt^ if) 
and 



\Gt,. (a {; f (•)))lli_«;[,T] < <T [t, T]) \\a / (■ 



l-a;[t,T] 



\a,oo;[t,T] 



Gi:^ if) 



< 



A.{g;[t,T]) ,2) 



a,\;[t,T] 



A 



l-2a 



<^ hi- J' 



\a,X;[t,T] 
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Now the inequalities (fTTI -i. ii) clearly follow, since for all A > 0, 

lk(-'/(-))lla,A;[i,T] 
-Ar 

Jt 



sup e 

re[t,T] 



< sup e 



Mo,T(l + |/(r)|) + Mo 



(r - 

•(r-sf +|/(r)-/(.) 



(r — s) 



l+a 



■ds 



< {Mo^T + Mo) 1 + 



^] 11+11^ ll-.A;[,T] 



Lemma 5 Let the assumption (Hi) be satisfied and < a < ^ A j3. Let f,h^ 



(t,T;M^) and g G W^-"''^{t,T;R''). Then fi 
follows 



or 



oo;[t,T] 



< R, 



oc;[t,T] 



Gt^ if) - (h) 



&fK{g-[t,T]) 



< 



a,A;[t,T] A^'^a 



< R it 



(12) 



/or a// A > 1, where 



A[t,T] (/) = sup 



(is. 



and C^^^ a constant only depending of Mq, Mj^, T, a, (3. 

Remark 6 //O < a < ||/|loo-[tT] — '^'^'^ ll/lli-a-[t t] — Gr , with Cr a constant 
independent oft, then 

A[i,T] (/) < Gn sup / \^ [ ds < , ^ , ,, Cn 



6-a{l + 6) 



re[t,T] Jt [r - s) 

Proof of Lemma [5l For the proof we use the ideas from [16j. By the inequality ([7]) 

we have 



Gt^ if) - GfJ (h) 



,A;[i,r] ^ A^-2° 
Remark that if |a;| , \y\ , \u\ , |f | < R, then 

2;) — a^'i (r, u) — 0"*'-' (s, y) + a^'^ (s, f 



< "'"^ Ik / (■)) -^i;h (■))IL,A;[,T] 



• (13) 



< 



[ {x-u,V^a''^{r,9x + {l-6)u)de- [ {y - v,V^a''^{s,9y + {1 - 9)v)) d9 
Jo Jo 



{x-y-u + v, V^a''^{s, 9y + {1 - 9)v)) d9 



+ 



(x - M, V^a^'J(r, + (1 - ^)m) - V^a''J(s, 9y + {I ~ 9)v)) d9 



and therefore 

\a (r, x) — cr (r, u) — a {s, y) + cr (s, v)\ < Mq \x — y — u + v\ 

+ Mr \x — u\ {^s — + \x — yf + \u — vf^ . 



Hence for \\f\\^^^t^T] < ^ and \\h\\^^^^^^^ < R 



\a{-J{-))~a{-,hm.Mt, 



sup e 



-Xr 



+ 



\(^irj{r)) -a{r,h{r))\ 

^ \a (r, / (r)) -a{r,h (r)) -a{sj (s)) + a{s,h {s))\ 



(r — s) 



l+a 



ds 



< (Mo + Mr) sup e-^'' 
+ |/(r)-/i(r)| 



/3 — a 



(r - 



where 



A[t,r] (/) = sup / 

■re[t,T] 



\f{r)-f{s)\' 



[r — s) 



vl+a 



and similar for h. We conclude that 

||a(-,/(-))-a(-,M-))IL,,,[,^] 



< (Mo + Mr) (^1 + ^— -J [1 + A[,,r] (/ + A[,,r] {h))] 11/ 
and the inequality ( fT2l) now follows from ( fT3l) . 



,A;[t,r] 



We also give similar estimates for 

fSU)= / &(^,/W)rf^- 



where h satisfies the assumptions (-f^2)- Very similar estimates are given in the paper 
of Nualart&Ra§canu [16] . 



Lemma 7 Let f e iy"'°°(t, T; M'^). Then F^^^ (/) = j^h{yj{y))dy e C^-" {t,T-W^) 
and for all X > 1 : 



(j) \\Fi!^{f) 

in) F,^!^ if) 



l-a;[t,T] 



<Ct^ fl + ll/lloo;[,T]), 



< 



(62) 



:,A;[t,T] A'= 

10 



Q,A;[t,T] 



(14) 



where Cq^^ and Cq'^'' are positive constants depending only on a,T and Lq 



Iff,he T; R'^) such that 



oo;[t,T] 



< R, 



oo;[t,T] 



< R, then 



Ff;^ if) - Fi:.^ (h) 



< 



c 



(b3) 
R 



Wf-hl 



[s — r) . 



,X;[t,T] - X- ■"'l".A;[i,T] 

for all A > 1, where C^^'^^ constants depending only on a,T and from (-^2)- 
Proof. It is easy to see that f/''^ (/) G ([t, T]) and for t < r < s < T 

i^ff (/) - FS (/)| = \FS (/)| < ^0 (1 + II. |loo;[,T], 

Hence the inequahty ffTH-j) follows with Cq^^^ = Lq (T + T*^ 
Denoting 

Ft,s if) = J J (r) dr 

we remark that 



(15) 



[s — r 



\a+l 



< \f{r)\dr + {s-r)-^-'{j;\f{u)\du)dr 

< f \f{r)\dr + -!\s-t)-- \f{u)\du 



and therefore 



Ft':' if) 



< T° + - sup e-^'{s-ty / \b{r,f{r))\dr 
»Mt,T] V "/ se[t,T] Jt 



<Lo(T" + -) sup [ e 

OiJ s&[t,T] Jt 



(s — r)° 



dr. 



1\ 1 T^-2a 
<Lo(T" + -l — . — (1 



a J A" 1 - 2a 



a,X;[t,T] 



I 1 A 

that clearly yields (IT^jj) with Cg^^^ = Lq ( H — ) since e 



a / 1 — 2a 



[s — r] 



-2a 



1 

A^ 

Let now f,he iy°'°^(t, T; M'^) and |/| < R and < R. Then as here above, for 
all A > 1, 



fI:^ if) - Fi:> {h) 



< (t° + -) sup e^^'is-ty f \b(rj(r))-b(r,h(r))\dr 
«:A;[i,r] V "/ se[t,T] Jt 



< r° + - 



1\ 1 T 



l-2a 



a A" 1 - 2a 



\\a,X;[t,T] ■ 



11 



Hence the inequality ([l5])holds with C^^^^ = ( T" + 



1\ T 



l-2a 



a J 1 — 2a 

m 

Finally we present some auxiliary estimates used in the sequel. 

Lemma 8 Let the assumptions (Hi) and {H2) be satisfied and < a < /3 A ^. IfYis 
a Holder continuous function with \\Y\\i_a-[tT] — ^ then there exist some positive con- 
stants = (R + l + T) Ln, dj^^ = (i?, Mq, T, a, p, A, (g)), ^ = 2 (1 + /?) Lq 
and C'^^ = C(4) (R^ Mq, T, a, p, A„ (g)) such that for all < t < s < T : 



[b{r,Yr)-b{t,Y,)]dr 



< Cg^ {s - tf-'^ and 



[air,Yr)-ait,Yt)]dg{r) 



< 4'^ is - 1) 



l+min{/3-«,l-2a} 



and for all < t < t < s < T : 

[b{r,Yr)-b{t,Yt)]dr 

(d) I [a{r,Yr)-ait,Yt)]dg{r) 



< ^ (s - r) and 



< Cf is - r 



Proof, (a) We have 



[b{r,Yr)-b{t,Yt)]dr 



<{s-t) sup \b{r,Yr)-b{t,Yt)\ 

re[t,s\ 

<{s-t) sup LR{XYr-Yt\ + \r-t\^) 

r&[t,s\ 

<{s-t) Lr sup I [i? + (1 + T)] (r - t) 

< 4^ (s - t)'-" 



l-a\ 



(b) By the assumptions (Hi) and Lemma we have 



\a{r,Yr)-a{e,Y0)\<Mo 
<Mo 



\r-9f + \Yr-Ye\ 

\r - ef + R\r - e\ 



and 



[air,Yr)-ait,Yt)]dgir) 



< 



' a{r,Yr) - a {t,Y,) - - 

(r - tr + 




t Jt 



a{r,Y,)-a{d,Ye 



■dddr 



(16) 



(17) 
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< A„ (g) Mo 



(r - tf-" + R{r-t) 



l-2a 



+ 




dr 



t Jt 



[r — 



-2a 



dOdr 



< cf [{s - tf-^^' + (s - tf~^- + {s- tf-"^"-' + {s 

< (^1^2) _ ^.|l+mm{/3-a,l-2a} 

(c) We have, by {H2 — ii) , 



[6(r,F,)-6(t,ri)]^^^ 



< (s-r) sup \bir,Yr)-bit,Yt)\ 

re[T,s] 

<Lo(s-r) (2 + |K,| + |y;|) 
<2{1 + R) Lo(s-r) 



'c?) Using (IH]) we deduce 
j [a{r,Yr)-a{t,Yt)]dg{r 

a{r,Yr)-a{t,Yt) 



< Aa (g) 



< (g) MoJ^ 



dr + 



r - r 




■or(r,r,)-^(^,>'e) 



{r-tf (r-t) 

+ -ti 



T 'J T 



r - r 



r - r 




< C 



(4) 
(4) 



1-Q I rp\ 



+ Aa{g)Mol I ((r-<-"-^ + i?(r 



,-2a 



dOdr 



a 



I \1— o I / a I / \2— 2a 

s-r + s-rr + s-r 



< - 



3 Fractional viability. Main result. 

Consider the stochastic differential equation driven by the fractional Brownian motion 
B with Hurst parameter l/2<if<l,P — a.s. cuGfi, 

/sMt psVt 
b{r,Xl'^)dr + a{r,X';'')dB^, se[t,T], (18) 

where 
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• B = (-B*)j;.xi ! -B*! i = ^,k, are independent fractional Brownian motions with 

Hurst parameter H, - < H < 1, and the integral with respect to i? is a pathwise 
Riemann-Stieltjes integral; 

• Xq = (Xq)^^^ is a d - dimensional random variable defined in a complete proba- 
bility space P); 

• 6 : [0,T] X W^^W^ and a : [0,T] x R'^''^ are continuous functions. 

Definition 9 Let fC = {K (t) : t G [0,T]} be a family of subsets ofW^. We shall say 
that K, is viable for the equation ( fl^) if, starting at any time t G [0,T] and from any 
point X E K {t) , at least one its solution X*'^ G K (s) for all s G [t, T] . 

Definition 10 The family K, is said to be invariant for the equation [W\] if, for any 
t G [0,T] and for any starting point x G K {t), all solutions {X^'^ : s G [t, T]} of the 
fractional stochastic differential equation fT8\} have the property 

X'/ G K (s) for all s G [t, T] . 

Remark that, in the case when the equation has a unique solution (which is the 
case for the equation (fTSi) under the assumptions (Hi) and {H2)), viability is equivalent 
with invariance. 

Assume that the mappings b and a from the equations (fT8|) are satisfying (Hi) and 

Definition 11 Let t G [0,T] and x e K {t) . Let \ < I - a < H. We say that the 
pair {b {t, x) , a (t, x)) is (1 — a) —fractional — contingent to K {t) in (t, x) if there 
exist random variable h = /i*'^' > 0, a stochastic process Q = Q*'^' -.Qx [t,t + h] ^ R'^ 
and for every R > such that \x\ < R there exist two random variables Hr,Hr > 
and a constant 7 = 7_r, G (0,1) which are independent of (t,Ji) (the constants Hji, 
Hr, 7_r depend only on R, Lr, Mq^t, Mq, Lq, T, a, (3, Aq, {B^^ ) such that for all 

s,T e [t,t + h] 

\Q{s)-Q{r)\<HR\s-T\'-'' and \Q {s)\ < Hr\s - t\'^^ 

satisfying 

x + {s-t)bit,x)+a {t, x) [5f - 5f ] +Q{s)eK (s) . 



The main result (existence result and characterization of the viability) of our paper 
is the following 
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Theorem 12 Let K, = {K (t) : t G [0,T]} he a family of nonempty closed subsets of 
Mf^. Assume that the maps b and a from the equations l[T^) are satisfying (Hi), {H2) 
with \ < H <l,l- H < (3, 6> 

Let 1 — H < a < at). Then the following assertions are equivalent: 

• K is viable for the fractional stochastic differential equation ( I7g|) . i.e. for all t G 
[0, T] and for all x e K (t) there exists a solution X*'^ {cu, ■) G C^"" ([t, T] ; M'^') 
of the equation 

/sVt psVt 
b{r, X';'')dr + a{r,Xl^''')dB^ , sG[t,T], a.s.uen, 

and X*'^ G (s) , for all s G [t, T] . 

• For all t G [0,T] and all x G K (t) , {b{t,x) ,(J {t,x)) is {1 — a) —fractional 

—contingent to K (t) in {t,x) . 



on 



4 Deterministic approach 

Let arbitrary fixed {t,x) G [0,T] x M'^. Consider the deterministic differential equation 

Xf = X + 6(r, Xf )rfr + J\ (r, X*^) dg (r) , s G [t, T] , (19) 

where g G W^-°^'^{t,T;R'') and the coefficients 6 : [0,T] x R'^ ^ R'^ and (x : [0,T] x 
R*^ — > M.'^^'^ are continuous functions satisfying the assumptions {Hi), {H2). Let a be 
arbitrary fixed such that 

f 1 5 

< a < ao = ^ 



^2" '1 + 6 

Nualart and Rascanu proved in [16] that if the assumptions (Hi) and {H2) are satis- 
fied then the equation (fT9|) has a unique solution which is {1 — a) — Holder continuous. 
In the following Lemma we shall proof that the Holder constant of this solution has 
the form Cq (1 + |x|) , with Cq a positive constant depending only on Mq^t, Mq, Lq T, 
a, j3. 

Lemma 13 Let the assumptions (Hi) and {H2) be satisfied. If X*'^ is a solution of 
the equation / fJ^) then X*'^ is {1 — a) -Holder continuous and 

ll^*1li-«;[,r]<^o(l + kl) 
where Cq is a constant depending only on Mq^t, Mq, Lq T , a, (3, (g) ■ 
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Proof. By Corollary H] and Lemma [7| we have for all A > 1 

||x*'^T 



la,A;[i,T] 



< \x\ + 



< \x\ + 



+ 



a,A;[t,r] 



(b2) 







1 + X'-^' 



^N + ^(l+ll^'1L.A;[.,T] 

for A = Aq > 1 sufficiently large, 



'■Mm 



+ 



A. ((7;[0,T])C^' 



{a2) 



Ai 



-2a 



\x 



t.x I 



:,A;[t,T] 



(a2) 



\ l-2a 
^0 



^2- 



(remark that the constant Cg^^"* + Aq, (g) C'^'^^ is independent of A). Then we have 

ll^'1L,A;[,T] <2(l + kl) and <2(l + |x|)e^o^. 

On the other hand, using the same lemmas we get 



(a2) 



ll-a;[t,T] 



< U\ + 



l-a;[t,T] 



l--a;[t,T] 



< Ixl + G 



(bl) 



\oo-\t,T]) + (^) '^o''^^ {} + l|-^*''''IL,oo;[t,T] 
<N+(cr+AQ(,)Cr) (l + ||^'1L.oo;,,T]) 

< |x| + (C^^ + A. {g) C\,"'^) (1 + 2e^"^ + 2 |x| e^«^) 
<Co(l + |x|) 

Hence X*'^ is (1 — a) —Holder continuous with the Holder constant Cq (1 + \x\) where 
Co depends only on Mo,r, Mq, Lq T, a, [3 and A^ {g) . The proof is now complete. 

■ 

Assume that the maps h and a from the equations (fT9l) are satisfying {Hi) and 

Definition 14 (Tangency property) Lei t G [0, T] anc? x E K [t) . We say that 
the pair {b {t, x) ,a {t,x)) is (1 — a) —fractional g— tangent to K (t) in (t, x) if there 
exist h = h^'^ > 0, and two functions U = U^'^ : [t, t + /i] — > W^, U (t) = 0, and 
V = 1/*'^ : [t,t + h] R'^'"', V (t) = 0, and for every R> such that \x\ < R there 
exist two constants Dr, Dr > independent of (t, fi) such that for all s,t E [t,t + h\ 



\U{t)-U{s)\<Dr\t-s 



and \V{t)-V{s)\<Dr\t-s 



min{/3, 1—a} 



(20) 
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and satisfying 

ps ps 

x + [b{t, x) + U (r)] dr + [a{t, x) + V (r)] dg (r) e K (s) , for all s E [t,t + h]. 

Definition 15 (Contingency property) Let t E [0,T] and x E K [t] . We say that 
the pair {b [t, x) ,a {t,x)) is (1 — a) —fractional g— contingent to K (t) in (t, x) if there 
exist h = h^'^ > 0, a function Q = Q*'^ : [t, t + /i] — > M.^ and for every R > such 
that \x\ < R there exist two constants Gr,Gr > independent of (t, li) and a constant 
7 = 7k £ (0, 1) also independent of (t,h) (the constants Gr,Gr, •jr depend only on 
R, Lr, Mq^t, Mq, Lq T, a, P and (g) ) such that for all s,t E [t,t + h\ 

|Q(r) -g(s)| < G/jIt-sI^"" and IQ < |s - t|^+^ 
and satisfying 

X + {s - t)b{t,x) + a (t, x) [g (s) - g (t)] + Q (s) E K (s) , for all s E [t,t + h] . 
We can now state the main result of the section. 



Theorem 16 Let /C = {K (t) : t E [0,T]} he a family of nonempty closed subsets of 
W^. Assume (Hi) and {H2) are satisfied and 



1 — H<a<ao = min 




Then the following assertions are equivalent: 

(j) /C is G^"""— viable for the fractional differential equation ( fTPj) . i.e. for any t E 
[0,T] and for any starting point x E K{t), there exists a solution X*'^ (■) E 
Gi-"([t,r]; M^) of the equation 

X'/ = x+ b{r, Xl'^)dr + / a(r, X*'")d(/ (r) , sE [t, T], (21) 

such that X*'^ E K {s) , for all s E [t, T]. 

(jj) For all t E [0, T] and allx E K (t) the pair {b {t, x) , a {t, x)) is (1 — a) —fractional 
g— tangent to K (t) in {t, x) . 

{Hi) For all t E [0, T] and allx E K (t) the pair {b (t, x) , a (t, x)) is (1 — a) —fractional 
g— contingent to K (t) in (t, x) . 
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Proof. Let < £ < 1. We denote by Cr, ... denote a generic positive 
constant independent of e and depending only on R, Lr, Mq^t, Md, Lq T, a, (3 and 

K (g) . 

(j) ^ (jj) : 

Let t e [0,T] and x e K (t) be arbitrary fixed and X*'^ G C^"" ([t,r] ; M^) a 

solution of the equation (12T|) such that X*'"^ G (s) , for all s G [t, T]. Let -Rq > such 
that \x\ < Rq. Then by Lemma [T^ 



ll-a;[i,r] 



<R = Coil + Ro) 



with Co depending only on Mq^t, Mq, Lq T, a, P and Aq (g) . Let /i = min {T — t, 1} 
Then 

/•S PS 

X'f = x+ b{r, X';'')dr + / a{r, X'/)dg (r) G K (s) , V s G [t, t + /i] . 
We clearly have for all s G [t, t + /i] 

/S PS 
[b{t,x) + U{r)]dr + [a{t, x) + V (r)] dg (r) 

where 

f/(r) = 6(r,X;'^) and V (r) = a(r, X*'^) - a (t, x) . 



Clearly U and F satisfy 
(jj) ^ (jjj) : 

Let I a; I < R. We verify that 

Q{s) = J^ U{r)dr + J^ V (r) dg (r) 

satisfies the Holder conditions from the definition of the contingency property. Indeed 
by we have 



U (r) dr 



and 



V{9)dgi9) 



< 



[U (r) - U (t)] dr 



[V (9) - V (t)] dg (9) 



< Dn \s-t\ 



2-a 



. f \ViO)-Vit)\ , nVi9)-Viu)\ 



't \ (o-t) 

< _ ^^l+mm{/3-a,l-2«} 
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Hence 



\Q {S)\ < {S - i)l+-M/3-«,l-2a} _ 



Now we prove that Q is (1 — a) —Holder continuous on [t, i + /i] . Let t < r < s < t+h. 
We have 



V{9)dg{9) 



< J [v{e)-v{t)]dg{e) 

< Kig) 



f \v{0)-v{t) 

V (0- tY 



\v{e) - v{u) 



du I d9 



and therefore 



\Q{s)-Q{r)\< 



U (r) dr 



V{9)dg{0) 



< C^' is - r) 



l-a 



(jjj) 



Let us fix t G [0, T], a; G (t) and < £ < 1. Let 7?o > be such that |a;| < Rq. 
We denote by As {t, x) the set of pairs (Tx, X) where Tx G [0, T] and X : [t, Tx] 
M.^ is a Holder continuous function satisfying 

(1) Xt = X, Xg ^ K (s) for all s G [t,Tx] , and there exists a positive constant 
-Bo > Rq depending only onit!o, L^^, Mo,t, M), Lq T, a, (3 and {g) , such that 

ll^lll-a;[t,Txl ^ ^0 

(2) The error function ^ : [t, Tx] 

^(s)^Xs-x- j^b(r, Xr) dr- j^a (r, Xr) dg (r) , s G [i, Tx] , 

satisfies 

(a) <£(s-i), for alls G 

(b) \C{t)-C{s)\<Do\t-s\'-'', foralls,TG [i,Tx], 

where the constant Dq depends only on Rq, Lji^, Mq^t, Mq, Lq T, a, j3 and Aq [g) . 
Remark that Bq and Dq are independent of e. 

The set As {t, x) is not empty because we can find {t, x) G As {t, x) . 
As (t, x) is an inductive set for the order relation 



{Tx,,X^{.))^{Tx,,X2{-)) 
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defined by 



Tx, < Tx, and (•) 



Zorn's Lemma impfies that there exists a maximal element {T'^,X'^) G Ac {t,x) . We 
shall prove by reductio ad absurdum that — T. 

Assume that < T. Denote X^s — x^. We have ||-'^^||i_Q-[t Tx] — -^o 
particular 

We know from the hypotheses that {b (T^, x^) , a (T^, x'^)) is (1 — a) —fractional contingent 
to K in {T'^^x^) , i.e. there exist /i^ > sufficiently small (for moment < /i^ < 
(T-T') A 1), :_[T^T^ + ^e] ^ K'', two constants Go = Gbo, Go = Gb,, > 
independent of (T^, /ig) and a constant 7 = e (0, 1) also independent of (T^, he) 
(the constants Gq,Gq, 7 depend only on Rq, Lfi^, Mq^t, M), -^0 o;, /3) such that for 
all s,r e [T^T^ + y 

\Q' (r) - (s)| < Go |r - and (s)| < Go |s - T^|'+^ 

and satisfying for all s e [T^,T^ + he] 

x' + {s- T') b (T^ x') + a (T^ x') [g (s) - g (T^)] + (s) E K (s) , 
We set — T^ + he and we define X^ : [t, S^] K as a, extension of X^ by 

r X-(s), ifse [t,T^], 

X" (s) = J + (s - T^) 6 (T^ x^) + a (T^ x') {g (s) - c/ (T^)) + g= (s) , 
I if s e [T^ 5^] . 



We will prove that the extension yS'^,X'^j G Ae {t,x). 

Step 1: Clearly Xf = x and X| e (s) for all s G [t, T] . 

Let us show that X^ 

Mo,T, Mo, Lq T, a, (3 and {g) . 
Let < s < r < S^. Then 



< i^Q^"* where Bq"^' > Rq and i^g'^'' depends only on Rq, L^^, 



(1) 



ATf - AT! 



< |t - s| \b (T% x^) I + Ac, (y) r (a) \a (T^ |t - + |g^ (r) - {s 



< |t-s|Lo(1+ + 



Aa (g) 



Mo,t{1 + \x'\)\t-s\'-" + Gr\t-s 



a 



\l-a 



< Gro \t-S 



1-a 



(22) 
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and 



Xi 



< 



x^ — X^E 



< C/j„Ti-° + i?o- Hence 



X' 


< 


X' 


+ 


X' 










l-a;[t,T^] 







Step 2: the error function. 

Let the error functions C,^ : [t, T 



^ and t ■■ [t, S'] R'^, 
r {s)=Xl-x- j^h{r,Xl)dr- j^a{r,Xl)dg{r). 
t {s)=Xl-x- j^h (r, xi) dr - j\ (r, X^) dg (r) . 



Clearly 



<e{s-t) for all s e [t,T']. 

Let s G [T^, 5*^] . Using Lemma [HI (the inequalities (fT6l) with Yr = X^ t = 
we have 





< 


— X — j 









x^ — X — b{r, X^) dr — a (r, X^) dg (r 



s -T')b{T',x') +a{T',x'){g{s) - g{T')) - / b[r,Xndr- / a ( r, ) rf^ (r) 



+ 

+ \QHs)\ 

<e{T' - 1) 

+ Go(s- T 



b{T',x') -b[r,Xf. 



rpe L 
£\l+7 



dr + 



0" (T^, x^) — 0" ( r, ) (r 



< £ (T^ - t) + eg,; (s - T'f-" + C'il {s - re^i+--{/3-.i-2"} + Go (s _ r^)i+T 

< e (T^ - t) + £ (s - T=) = £ (s - t) 
for hs sufficiently small such that 

Hence 



<e{s-t) for all s G [t, S' 



Let now < r < s < S"^. Then by Lemma [8] (the inequalities (fT7|l with y,. = X, 

t = T'^ and = X|.e we have 



e (^) - e (r) 



< 



6(T^a;^) -6(r,X,' 



(ir + 



(T(T^x^) - a (r,X," 
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< C^S - ^) + - ^)'"" + Go\r- s\ 



vl-a 



<Cn,{s-T) 
From the definition of Ae {t, x) , for r, s G [t, T*^ 



r («) - r (r) = If (s) - f (r)| < L'o - r)^"" , V r, s e [t, . 
We conclude 

r («) - r (r) I < (C«„ V L>o) |r - , V r, s e [t, S^] . 

We arrived to prove that ^5"^, X^^ is proper extension of {T^, X^) , that contradicts 

the maximaUity of (^^ X^) in (t, x) . Therefore = T. 

Let (T, X^) be a maximal element of Ae {t,x) . Then from the definition of Ae {t, x) 
we have Xf = x, X| G (s) for all s G [t, T] , and there exists a positive constant 
-Bo > -Ro depending only on Rq, Lr^, Mq^t, Mq, Lq T, a, j3 and {g) , such that 



The error function : [t, T] 



e{s)=XI-x- / b{r,X'^)dr- / a {r, X'^) dg {r) 
Jt Jt 



(23) 



satisfies 



(a) ir (s) I < £ (s - t) , for all s G [t, T] 



(b) \e (r) - {s)\ <Do\t- , for all r G [t, T] , 

where the constant Dq depends only on Rq, Lr^, Mq^t, Mq, Lq T, a, (5 and {g) 
We now estimate oo [t r] ■ A > 0. Then 

lll^^llla,A;[t,r] — ll^^lla,oo;[t,T] 



= sup + 



se[t,T\ L Jt \s- r) 



-dr 



<e(T-t)+ sup Tie-ei^ 



o: |Ss Sr I 



s — r 
s ni+" 



-dr 



<e{T-t) + [2e (T - t)]^"^ sup 



s — r 



(s — r) 



l+a 



dr 



^e{T-t) + [2e (T - t)]^" (T - 
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It remains now to prove that the hmit of the sequence exists as £ ^ and this 

hmit is a solution to the differential equation f|T9l) 

Let Q < e,rj <\. Using the estimates ( |T2|) and ( |T5|) . we get 



< 



< 



X[t,T] 



+ iiie-e'iiL,A;[,T] 



gr i^^ - ^1IL.A;[,Ti+ (1 + A[t,T] (x^) + A[,^] (x^)) iiix^ - x^iiL 



X;[t,T] 



Since ||X''||^_^.jj < Bq , < a < then 

A[i,T] (X^) = SUp^,g[i_r] / — <BoSUp^^[tT] 

Jt (r - s) Jt 



[r — s 



■a)5 



r — S 



\a+l 



-ds 



< 



(5 -a (1 + 5) 
5-a{l + 5) 



and therefore 



111^^ - X1IL,A;[*, 



< 



4'^ 

Ho 



+ + Cr, r/^- 



«,A;[t,r] 



Let A = A > 1 such that 



Ho 



6-a{l + 6) VA" Ai 



-2a 



We deduce 



1 

< -. 



11^^ - ^''lloo;[,T] < e'^ |||X^ - < 2CR,e'^ (^i"" + 



Hence there exists X*-^ such that X= X*-^' in C([t,T];/s:) and X^ ^ X*'^ in 
T; M^) as £ ^ 0. Since for all s, r G [t, T] : 



||X^ 



loo;[t,T] 



^f-xy 

+ : < i^o 
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then passing to limit as e — we obtain 

II Vt,x II ^ D 

F \\l-a;[t,T] ^ ^0- 

Since X| G K (s) for all s G [t, T] , clearly follows, as e — > 0, that 

XsEK (s) for all s G [t, T] . 

Passing to limit in fl23l) we infer that X*'^' is a solution of the differential equation f[T9l) 
starting at t from x and evolving in the tube {{s,y) : s E [t,T] , y E K (s)} . 
The proof is complete. 



5 Proof of the main result 

Let be fixed a parameter 1/2 < H < 1. Consider B = {Bt,t G [0,T]} be a M.^ valued 
fractional Brownian motion with parameter H defined in a complete probability space 
(fi, J^, P) . From ([IDit follows that 

E{\Bt-B,f) = It-sf"" 
and as a consequence, for any p > I, 

\\B, - B^l = (E {\B, - Bf)f^ = cp\t-sf. 
It is known that the random variable 

g= P,/ , sup \{Dlz''Br-)is)\ 

T [1 - a) t<s<r<T 

has moments of all order. As a consequence. As a consequence, if u = {uf, t G [0,T]} 
is a stochastic process whose trajectories belong to the space W^"'^(0, T; M'^^'^), with 
1 — H < a < the pathwise integral UgdBs exists in the sense of Definition [2] and 
we have the estimate 



T 

U.4B, 







<G \\u\\,. 



Moreover, if the trajectories of the process u belong to the space iy"'°°(0, T; M'^^'^'), 
then the indefinite integral Ut = j^UsdBg is Holder continuous of order 1 — a, and 
the estimates from Proposition [3] hold. 

Proof of the Theorem 1121 Considering the previous observations, the solution 
follows directly from the deterministic Theorem | 
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